Observation of Excited States in a Graphene Quantum Dot 



S. Schnez, F. Molitor, C. Stampfer, J. Giittinger, I. Shorubalko, T. Ihn, and K. Ensslin 

Solid State Physics Laboratory, ETH Zurich, 8093 Zurich, Switzerland 

We demonstrate that excited states in single-layer graphene quantum dots can be detected via 
direct transport experiments. Coulomb diamond measurements show distinct features of sequential 
tunneling through an excited state. Moreover, the onset of inelastic co-tunneling in the diamond 
region could be detected. For low magnetic fields, the position of the single-particle energy levels 
fluctuate on the scale of a flux quantum penetrating the dot area. For higher magnetic fields, the 
transition to the formation of Landau levels is observed. Estimates based on the linear energy- 
momentum relation of graphene give carrier numbers of the order of 10 for our device. 



Graphene [H, 0], the first real two-dimensional (2D) 
solid, consists of a hexagonal lattice of carbon atoms 
providing highly mobile electrons [3, 0] for future ap- 
plications in electronics, spintronics 5[ and information 
processing [6]. However, confinement of charge carriers 
in graphene cannot be achieved as easily as in conven- 
tional two-dimensional electron gases by using electro- 
static gates because of the gapless nature of graphene 
[1] and a relativistic phenomenon called Klein tunnel- 
ing 0, H[ . Cutting graphene into a desired geometry is 
an alternative to overcome this obstacle. Well-controlled 
nanostructures, such as nanoribbons quan- 
tum interference devices [12, UjS, Lb! , and single-electron 
transistors [H, EE 1 1 71 ] have been created in several labs 
to date. Small spin-orbit and hyperfine interactions have 
been theoretically predicted [l8[, promising spin deco- 
herence times superior to the GaAs material system in 
which solid-state spin qubits are most advanced today 
[H, Eq[. Therefore, the identification of individual or- 
bital quantum states, well established in GaAs quan- 
tum dot devices, have so far remained on the wish list 
of physicists aiming at quantum information processing 
with graphene. 

An atomic force microscope image of our quantum dot 
(QD) is shown in Fig. la. It was fabricated with the 
standard procedure: Mechanical exfoliation of natural 
graphite led to single-layer graphene flakes. The desired 
structure was defined with electron beam lithography and 
subsequently cut using reactive ion etching based on Ar 
and O2. Contacts were also defined with electron beam 
lithography; then gold contacts were evaporated on top 
fTHd . The single layer quality was experimentally verified 
with Raman spectroscopy 21]. The QD device consists 
of two about 60 nm and 70 nm wide graphene constric- 
tions separating source (S) and drain (D) contacts from 
the graphene island (diameter 140 nm). The island can 
be tuned by a nearby plunger gate (PG), whereas the 
overall Fermi level is adjusted with a highly doped sili- 
con back gate (BG). The sample was annealed for about 
24 hours at 400 K directly before cool down. The experi- 
ments were carried out in a dilution refrigerator at a base 
temperature of 40 mK. Measuring the current I through 
the QD as a function of back gate voltage Vb g allows us 
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FIG. 1: (Color online) (a) Atomic force micrographs of the 
measured quantum dot. The quantum dot can be tuned by a 
nearby plunger gate (PG). The central island is connected to 
source (S) and drain (D) contacts by two constrictions. The 
diameter of the dot is 140 nm. (b) A backgate sweep shows a 
transport gap from roughly Vbg = — 8 V to 8 V (bias voltage 
of Vbias = 3mV). (c) Coulomb blockade measured with the 
back gate at an electronic temperature of 200 mK and a bias 
voltage of Vbias = 16 /xV. 



to identify a transport gap [15[ extending roughly from 
Vbg = — 8V to 8V (Fig. lb). Since the gap is centered 
around zero back gate voltage, we have little doping of 
our graphene device. Characteristic peaks in the gap- 
region were identified as Coulomb peaks (Fig. lc) and 
could be used to estimate an upper bound for the elec- 
tronic temperature. This was found to be around 200 
mK. In the following measurements, we set the back gate 
voltage to zero in order to tune the device close to the 
charge neutrality point. 

Coulomb diamond measurements [22[, i.e., plots of the 
differential conductance G — d//dVbi as , as a function of 
the quantum dot bias voltage Vbi as and plunger gate volt- 
age Vp g are shown in Fig. 2. Within this plunger gate 
voltage range, no charge rearrangements were observed 
and the sample was stable for more than two weeks. We 
extract a typical energy scale of the order of 10 meV. A 
strong fluctuation of the addition energy over the plunger 
gate voltage range —0.1 V < V pg < 1.2 V (full data range 
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FIG. 2: (Color online) (a) Differential conductance G (loga- 
rithmic plot) as a function of source-drain voltage Vbias and 
plunger gate voltage V pg . (b) The right panel is a zoom of the 
enframed area in (a) . An excited state is clearly visible (white 
arrow). A cut along the dashed line at Vbias = — 2.87mV 
is shown in the left panel (here G is measured in units of 
10~ 3 e 2 /h and was smoothed over 4 points), (c) Stability dia- 
gram at different plunger gate voltages. Several excited states 
are visible as shown by dashed lines. In the upper diamond, 
regions of higher conductance (in the left part of the diamond) 
can be seen. This is interpreted as a signature of co-tunneling 
in a graphene quantum dot (see arrow). In all measurements 
shown in Figs. 2 and 3, the back gate voltage was set to 
Vb g = V and the electronic temperature was around 200 
mK as deduced from the Coulomb peak width. 



not shown), corresponding to an energy range of around 
100 meV is observed, indicating the importance of quan- 
tum confinement effects. 

This is supported by the observation of excited states, 
which appear in Fig. 2a as distinct lines of increased 
conductance running parallel to the edge of the Coulomb 
diamonds [22[. Fig. 2b showing a close up of Fig. 2a 
allows to extract an excitation energy A « 1.6 meV as 
marked by the white arrow. A line cut at Vbi as = 2.78 mV 
(dashed line) presented in the left panel of Fig. 2b 
shows the peak of the excited state at finite bias (ar- 
row). The broadening of the peak significantly exceeds 



thermal broadening and might be due to the energy de- 
pendent coupling of the excited state to the graphene 
leads. 

Fig. 2c shows two Coulomb diamonds at lower plunger 
gate voltage, where more than one excited state is ob- 
served as a function of increasing energy, as shown by 
pairs of dashed lines. These excitations are found at ener- 
gies of around 1.6 meV and 3.3 meV (black dashed lines) 
and 2.1 meV and 4.2 meV (white dashed lines), respec- 
tively. The observation of excitations at finite source- 
drain voltage finds support by the detection of inelas- 
tic co-tunneling onsets at lower bias. Inside the upper 
Coulomb diamond of Fig. 2c, we distinguish between 
regions of suppressed and slightly elevated conductance 
separated by the dotted line. The edge of this conduc- 
tance step is aligned with the (first) excited state outside 
the diamond at an energy of 1.6 meV as highlighted by 
an arrow. 

The number of charge carriers on the quantum dot 
can be estimated by using the linear density of states of 
graphene D(E) = 2E / (jrvph 2 ) , where vp ~ 10 6 m/s is 
the Fermi velocity [2[ . We find for the single-particle level 
spacing A(7V) = hvp/^d^/N) [26], where N is the num- 
ber of charge carriers and d is an effective dot diameter 
related to the dot area A according to A = ird 2 /4. From 
the measured excitation energy and the lithographic di- 
mensions of the dot (d = 140 nm) we estimate the num- 
ber of charge carriers on the dot to be of the order of 
10. This estimation obviously changes if the dot area 
is larger, e. g. if the dot was not defined by the two 
constrictions. The number of charge carriers would then 
be even lower. However, magnetic field sweeps corrob- 
orate our estimation of the dot size as pointed out be- 
low. The current was below measurement resolution at 
Vp g < —0.5V so that smaller charge carrier numbers and 
the potential electron-hole crossover could not be stud- 
ied. The Coulomb diamonds shown in Fig. 2c scatter 
significantly stronger in size than those presented in Fig. 
2a. This might be a consequence of the lower number of 
charge carriers on the quantum dot. According to theo- 
retical considerations [23], electron- hole symmetry could 
manifest itself by an enhanced confinement energy. This 
should be detectable by the type of experiment described 
here. 

Using typical values for the addition energy and ex- 
citation energies we estimate the charging energy to be 
~ 8.5 me V, which agrees reasonably well with the en- 
ergy estimated from a disk model AEc = e 2 /(4ee e ff<i) « 
12 meV. We assumed the effective dielectric constant in- 
cluding vacuum and the Si02 to be e e ff = (l + 4)/2 = 2.5. 

In order to further explore the excitation spectrum, 
we show the energy shift of nine consecutive Coulomb 
peaks in a magnetic field applied normal to the plane 
of the quantum dot in Fig. 3. The vertical energy axis 
was obtained by converting plunger gate voltage into en- 
ergy usign the measured lever arm (a pg = 0.075). In 
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FIG. 3: (Color online) Experimental energy spectrum of the 
quantum dot in a perpendicular magnetic field. The typical 
magnetic field scale at which a significant change is expected 
is approximately given by one flux quantum <E>o = h/e per 
dot area, i.e. 4$o /nd 2 = 270 mT and is indicated by the 
black arrow. Starting around B = 4 T, a regular pattern with 
characteristic linear slopes evolves (see dashed lines) which 
shows the transition from single-particle fluctuations to B- 
field dependence. 



the constant-interaction model, the ground-state energy 
of an TV-particle quantum dot can be written as the sum 
of the single-particle energies £i(B) plus an electrostatic 
charging energy NEc- The ground-state energy is tuned 
by the gate voltage V g . The experiment was done in the 
zero-bias regime; hence we measured the chemical po- 
tential of the TVth Coulomb resonance as explained in 
[23( | . Experimentally, the single-particle energy Sn(B) 
of the Nth Coulomb resonance is then determined by 
e N (B) = ea pg V g res (N,B) + NE C + const. The con- 
stant part and the electrostatic contribution NEc are 
subtracted such that consecutive peaks touch each other 
( alter natingly shown as red triangles and blue circles, 
respectively). Characteristic lines (see dashed lines in 
Fig. 3) linear in B with slopes of around ±2.5meV/T 
can be seen. This strong 5-field dependence cannot be 
explained by the Zeeman effect, which would result in 
a slope of gfiB = 116/ieV/T, assuming a ^-factor of 
g = 2. For higher magnetic fields, the Landau level 
degeneracy increases and fewer Landau levels are filled. 
Consequently, the energy spectrum is expected to evolve 
from single-level fluctuations into a regular pattern. This 
transition can be seen at around 4 T. Recent theoretical 
calculations are in qualitative agreement with our exper- 
imental data [23|, l24| . 

Quantum dots are envisioned as possible building 
blocks for a future quantum information processor [25[. 
The preparation and detection of well-defined orbital and 
spin states are an essential prerequisite for this purpose. 
On the other hand, graphene might be a highly suitable 
material for spin manipulations in a condensed-matter 



environment because of its expected long spin-coherence 
times In this paper, we showed that the direct mea- 
surement of excited states in graphene quantum dots 
through transport experiments is possible. This is a 
first and essential step towards possible experiments with 
graphene quantum structures. 

Acknowledgement: We thank K. S. Novoselov and 
A. K. Geim for valuable discussion. Financial support by 
ETH Zurich and the Swiss Science Foundation is grate- 
fully acknowledged. 



[1] A. K. Geim and K. S. Novoselov, Nat. Mat. 6, 183 (2007). 
[2] A. H. Castro Neto, F. Guinea, N. M. Peres, and A. K. 

Geim, Rev. Mod. Phys. (in press), cond-mat 0709.1163vl 

(2007). 

[3] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 
A. A. Firsov, Nature (London) 438, 197 (2005). 

[4] Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 
(London) 438, 201 (2005). 

[5] N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkmann, 
and B. J. van Wees, Nature (London) 448, 571 (2007). 

[6] B. Trauzettel, D. V. Bulaev, D. Loss, G. Burkard, Nature 
Phys. 3, 192 (2007). 

[7] N. Dombey and A. Calogeracos, Phys. Rep. 315, 41 
(1999). 

[8] M. I. Katsnelson, K. S. Novoselov, A. K. Geim, Nature 

Phys. 2, 620 (2006). 
[9] Z. Chen, Y. Lin, M. Rooks, and P. Avouris, Physica E 

40, 228 (2007). 

[10] M. Y. Han, B. Ozyilmaz, Y. Zhang, and P. Kim, Phys. 
Rev. Lett. 98, 206805 (2007). 

[11] X. Li, X. Wang, L. Zhang, S. Lee, and H. Dai, Science 
319, 1229 (2008). 

[12] F. Miao, S. Wijeratne, Y. Zhang, U. C. Coskun, W. Bao, 
and C. N. Lau, Science 317, 1530 (2007). 

[13] S. Russo, J. B. Oostinga, D. Wehenkel, H. B. Heersche, 
S. S. Sobhani, L. M. K. Vandersypen, and A. F. Mor- 
purgo, Phys. Rev. B 77, 085413 (2008). 

[14] F. Molitor, L. Hiifner, A. Jacobsen, A. Pioda, K. Ensslin, 
and T. Ihn, unpublished. 

[15] C. Stampfer, J. Giittinger, F. Molitor, D. Graf, T. Ihn, 
and K. Ensslin, Appl. Phys. Lett. 92, 012102 (2008). 

[16] L. A. Ponomarenko, F. Schedin, M. I. Katsnelson, R. 
Yang, E. H. Hill, K. S. Novoselov, and A. K. Geim, Sci- 
ence 320, 356 (2008). 

[17] C. Stampfer, E. Schurtenberger, F. Molitor, J. Giittinger, 
T. Ihn, and K. Ensslin, Nano Letters 8, 2378 (2008). 

[18] H. Min, J. E. Hill, N. A. Sinitsyn, B. R. Sahu, L. Klein- 
man, and A. H. MacDonald, Phys. Rev. B 74, 165310 
(2006). 

[19] J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, 
A. Yacoby, M. D. Lukin, CM. Marcus, M. P. Hanson, 
and A. C. Gossard, Science 309, 2180 (2005). 

[20] F. H. L. Koppens, C. Buizert, K. J. Tielrooij, I. T. Vink, 
K. C. Nowack, T. Meunier, L. P. Kouwenhoven, and 
L. M. K. Vandersypen, Nature (London) 442, 766 (2006). 

[21] D. Graf, F. Molitor, K. Ensslin, C. Stampfer, A. Jungen, 
C. Hierold, and L. Wirtz, Nano Letters 7, 238 (2007). 



4 



[22] L. P. Kouwenhoven, C. M. Marcus, P. L. McEuen, S. 
Tarucha, R. M. Westervelt, and N. S. Wingreen, Proceed- 
ings of the NATO Advanced Study Institute on Meso- 
scopic Electron Transport E345, 105 (1997). 

[23] S. Schnez, K. Ensslin, M. Sigrist, and T. Ihn, Phys. Rev. 
B 78, 195427 (2008). 

[24] P. Recher, J. Nilsson, G. Burkard, and B. Trauzettel, 



larXiv:0810.04l9l (unpublished) . 
[25] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

[26] This estimation does not hold for very small in partic- 
ular not for N = 0. However, it is a reasonable estimation 
for electron numbers exceeding N > 5. 



